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ABSTRACT 

We present sensitive 7 mm observations of the H53q; recombination line and 
adjacent continuum, made toward the Orion BN/KL region. In the continuum 
we detect the BN object, the radio source I (GMR I) and the radio counterpart 
of the infrared source n (Orion-n). Comparing with observations made at simi- 
lar angular resolutions but lower frequency, we discuss the spectral indices and 
angular sizes of these sources. In the H53a line we only detect the BN object. 
This is the first time that radio recombination lines are detected from this source. 
The LSR radial velocity of BN from the H53q; line, vlsr = 20.1 ± 2.1 km s~^, 
is consistent with that found from previous studies in near-infared lines. While 
the continuum emission is expected to have considerable optical depth at 7 mm. 
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the observed H53a line emission is consistent with an optically-thin nature and 
we discuss possible explanations for this apparent discrepancy. There is evidence 
of a velocity gradient, with the NE part of BN being redshifted by ~ 10 km s~^ 
with respect to the SW part. This is consistent with the suggestion of Jiang et 
al. that BN may be driving an ionized outflow along that direction. 

Subject headings: ISM: individual (Orion) — radio continuum: ISM — radio 
hues: ISM 



1. Introduction 

Located at about 1' to the NW of the Orion Trapezium, the BN/KL region has been, 
as the closest region of massive star formation, the subject of extensive studies. Recently, 
Rodriguez et al. (2005) and Gomez et al. (2005) reported large proper motions (equivalent 
to velocities of the order of a few tens of km s~^) for the radio sources associated with the 
infrared sources BN and n, as well as for the radio source I. All three objects are located at 
the core of the BN/KL region and appear to be moving away from a common point where 
they must all have been located about 500 years ago. Even when these proper motions 
are now available, there is no radial velocity information for these three sources, with the 
exception of the near-infrared spectroscopic study of BN made by Scoville et al. (1983), that 
report an LSR radial velocity of -|-21 km s~^ for this source. In this paper we present 7 
mm continuum and H53q; radio recombination line observations of the BN/KL region in an 
attempt to obtain additional information on the radial velocities of these sources. 



2. Observations 

The 7 mm observations were made in the B configuration of the VLA of the NRAoQ, 
during 2007 December 14. The central rest frequency observed was that of the H53a line, 
42951.97 MHz, and we integrated on-source for a total of approximately 3 hours. We ob- 
served in the spectral line mode, with 15 channels of 1.56 MHz each (10.9 km s~^) and 
both circular polarizations. The bandpass calibrator was 0319+415. A continuum channel 
recorded the central 75% of the full spectral window. The absolute amplitude calibra- 
tor was 1331+305 (with an adopted flux density of 1.47 Jy) and the phase calibrator was 
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0541—056 (with a bootstrapped flux density of 1.78±0.08 Jy). The phase noise rms was 
about 30°, indicating good weather conditions. The phase center of these observations was 
at q;(2000) = 05^* 35"^ 14fl3; 5(2000) = -05° 22' 26'.'6. 

The data were acquired and reduced using the recommended VLA procedures for high 
frequency data, including the fast-switching mode with a cycle of 120 seconds. Clean maps 
were obtained using the task IMAGR of AIPS with the ROBUST parameter set to 0. 

3. Continuum Analysis 

3.1. Spectral Indices 

In Figure 1 we show the image obtained from the continuum channel. Three sources, 
BN, I and n, are evident in the image. No other sources were detected above a 5-0" lower limit 
of 1.75 mJy in our 1' field of view. The positions, flux densities, and deconvolved angular 
sizes of these sources are given in Table 1. The continuum flux density of the sources has 
been obtained from the hne-free channels. The line emission will be discussed below. The 
flux density obtained at 7 mm by us for BN is in good agreement with the values previously 
reported in the literature: we obtain a flux density of 28.6±0.6 mJy, while values of 31±5 and 
28.0±0.6 were obtained by Menten & Reid (1995) and Chandler & Wood (1997), respectively. 
In the case of source I, the agreement is acceptable, since wc obtain a flux density of 14.5±0.7 
mJy, while values of 13±2 and 10.8±0.6 mJy were reported by Mcntcn & Reid (1995) and 
Chandler & Wood (1997), respectively. Careful monitoring would be required to test if the 
radio continuum from source I is variable in time. 

The spectral indices determined from our 7 mm observations and the 3.6 cm observations 
of Gomez et al. (2008) are given in the last column of Table 2. Our spectral indices for 
BN and I are in excellent agreement in this spectral range with the more detailed analysis 
presented by Plambeck et al. (1995) and Beuther et al. (2004). 

We have detected source n for the flrst time at 7 mm and this detection allows the 
flrst estimate of the spectral index of this source over a wide frequency range. The value of 
0.2±0.1 suggests marginally thick free- free emission, as expected in an ionized outflow. This 
supports the interpretation of this source as an ionized outflow by Gomez et al. (2008). The 
position given by us in Table 1 is consistent with the extrapolation of the proper motions of 
this source discussed by Gomez et al. (2008). 
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3.2. Deconvolved Angular Sizes 

The radio source I has parameters consistent with an optically thick free-free source 
(spectral index of 1.5 ± 0.1). Beuther et al. (2004) suggest that this spectral index is either 
the result of optically thick free-free plus dust emission, or H~ free-free emission that gives 
rise to a power- law spectrum with an index of ~1.6. 

In the case of the radio source associated with the infrared source n we only have 
an upper limit to its size at 7 mm. In addition, Gomez et al. (2008) report important 
morphological variations over time in this source that suggest that comparisons at different 
frequencies should be made only from simultaneous observations. 

In the case of BN, the frequency dependences of flux density and angular size (this last 
parameter taken to be the geometric mean of the major and minor axes reported in Tables 
1 and 2) can be accounted for with a simple model of a sphere of ionized gas in which the 
electron density decreases as a power-law function of radius, Ue oc r~°. In this case, the flux 
density of the source is expected to go with frequency d& Sj, oc i/{6-2-4Q;)/(i-2a) ^j^g angular 
size is expected to go with frequency as 6i, oc (Reynolds 1986). The frequency 

dependences of flux density {Sy oc t/i-i±0-i) and angular size [Ou oc //-0-36±o.i2^ ^^j, -qt^ 
consistent with a steeply declining electron density distribution with power law index of 
a = 3.0 ±0.3. The continuum spectrum of BN produced by Plambeck et al. (1995) indicates 
that a constant spectral index extends from 5 to 100 GHz. 



4. Analysis of the H53q; Recombination Line Emission 

4.1. Radial LSR Velocity 

We clearly detected the II53q; line emission only from BN. The spectrum is shown in 
Figure 2. The parameters of the Gaussian least squares fit to the profile are given in Table 
3. We note that the radial LSR velocity determined by us, -1-20.1 ± 2.1 km s"^, agrees well 
with the value of -|-21 km s~^ reported by Scoville et al. (1983) from near-IR spectroscopy. 
In a single dish study of the H41a line made with an angular resolution of 24" toward Orion 
IRc2, Jaffe & Martin-Pintado (1999) report emission with vlsr = -3.6 km s~^. Most likely, 
this is emission from the ambient H II region, since its radial velocity practically coincides 
with the value determined for the large H II region (Orion A) ionized by the Trapezium stars 
(e. g. Peimbert et al. 1988). The single dish observations of the H51q; emission of Hasegawa 
& Akabane (1984), made with an angular resolution of 33", most probably come also from 
the ambient ionized gas and not from BN. 
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4.2. LTE Interpretation 

If we assume that the hne emission is optically thin and in LTE, the electron tempera- 
ture, Tg*, is given by (Mezger & Hoglund 1967; Gordon 1969; Quireza et al. 2006): 




where vl is the line frequency, Sq is the continuum flux density, Sl is the peak line flux 
density, Af is the FWHM line width, and y'^ is the ionized helium to ionized hydrogen 
abundance ratio. In the case of BN, we can adopt y+ ~ given that the source is not of 
very high luminosity, and using the values given in Tables 1 and 3, we obtain T* ~ 8, 200 K. 
This value is similar to that determined for the nearby Orion A from radio recombination 
lines (e. g. Lichten, Rodriguez, & Chaisson 1979). 

It is somewhat surprising that we get a very reasonable estimate for T* when our 
previous discussion seemed to imply that BN is partially optically thick at 7 mm. One 
possibility is that we have two effects fortuitously canceling each other. For example, the 
optical thickness of the source will diminish the line emission, while maser effects (such as 
those observed in MWC 349; Martm-Pintado et al. 1989) will amphfy the hne. However, in 
an attempt to understand this result in LTE conditions, we will discuss the expected LTE 
radio recombination line emission from a sphere of ionized gas in which the electron density 
decreases as a power-law function of radius, rie oc r^". As noted before, the modeling of the 
continuum emission from such a source was presented in detail by Panagia & Felli (1975) 
and Reynolds (1986). The radio recombination line emission for the case a = 2 has been 
discussed by Altenhoff, Strittmatter, & Wendker (1981) and Rodriguez (1982). Here we 
generalize the derivation of the recombination line emission to the case of a > 1.5. This 
lower limit is adopted to avoid the total emission from the source to diverge. 

For a sphere of ionized gas, the free-free continuum emission will be given by (Panagia 
& FeUi 1975): 



Sc^2n^B,j^ (1 - expi-Tcm) C dC, (2) 

where tq is a reference radius, d is the distance to the source, B^, is Planck's function, ^ 
is the projected radius in units of ro, and tc{C) is the continuum optical depth along the 
line of sight with projected radius ^. On the other hand, the free-free continuum plus radio 
recombination line emission will be given by an equation similar to eqn. (2), but with the 
continuum opacity substituted by the continuum plus line opacity (Rodriguez 1982): 
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Sl+c = 27r^S. (l - expi-TL+cmyd^, (3) 

where tl^c{0 P^^^ continuum optical depth along the line of sight with projected 

radius ^. 



The line-to-continuum ratio will be given by: 

Sc Sc 



(4) 



The opacity of these emission processes depends on projected radius as (Panagia & FeUi 
1975): 

riO^r^"'-'^. (5) 
We now introduce the definite integral (Gradshteyn & Ryzhik 1994) 



f°° 1 2 2 

/ [1 — exp(— iJ,x~^)] X dx — — jip r( — ), 
Jo P P 



(6) 



valid for > and p > and with F being the Gamma function. Substituting eqns. (2) 
and (3) in eqn. (4), and using the integral defined in eqn. (7), it can be shown that 



Sl 
Sc 



Kc 



l/(a-0.5) 

- 1, (7) 



where kl and kc are the line and continuum absorption coefficients at the frequency of 
observation, respectively. In this last step we have also assumed that the opacity of the line 
and continuum processes are proportional to the line and continuum absorption coefficients, 
respectively, that is, that the physical depths producing the line and continuum emissions 
are the same. Under the LTE assumption, we have that 



Kc \GHzJ \kJ \kms-^J ^ ^ ' ^' 



For 1/ < 43 GHz and typical parameters of an H II region, we can see from eqn. (8) 
that kl < Kc, and eqn. (7) can be approximated by: 
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Sl 1 {k-l 

Sc (a - 0.5) [kc. 



(9) 



That is, the expected optically-thin, LTE hne-to-continuum ratio: 



Sc 



(10) 



becomes attenuated by a factor l/{a — 0.5). In the case of a = 2, the factor is 2/3, and we 
reproduce the result of Altenhoff, Strittmatter, & Wendker (1981) and Rodriguez (1982). In 
the case of BN, we have that a ~ 3, and we expect the attenuation factor to be 2/5. If BN 
can be modeled this way, we would have expected to derive electron temperatures under the 
LTE assumption (see eqn. 1) of order 



T*{a = 3) ~ 2.2 T*{thin). (11) 

However, from the discussion in the first paragraph of this section observationally we 
determine that 



T*{a^3) c^T*{thin). (12) 

Summarizing: i) BN seems to have significant optical depth in the continuum at 7 mm, 
ii) this significant optical depth should attenuate the observed recombination line emission 
with respect to the optically-thin case, but iii) the line emission seems to be as strong as in 
the optically-thin case. 

As possible explanations for the "normal" (apparently optically-thin and in LTE) radio 
recombination line emission observed from BN we can think of two options. The first is 
that, as noted before, there is a non-LTE line-amplifying mechanism that approximately 
compensates for the optical depth attenuation. The second possibihty is that the free-free 
emission from BN at 7 mm is already optically thin. However, this last possibility seems to 
be in contradiction with the results of Plambcck ct al. (1995) that suggest a single spectral 
index from 5 to 100 GHz. Observations of radio recombination lines around 100 GHz are 
needed to solve this problem. 

A comparison with the H53a emission from the hypercompact H II region G28.20-0.04N 
is also of interest. The continuum flux densities from this source at 21, 6, 3.6, and 2 cm are 
49, 135, 297, and 543 mJy, respectively (Sewilo et al. 2004). At 7 mm the continuum flux 
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density is 641 mJy (Sewilo et al. 2008), indicating that the source has become optically thin 
at this wavelength. Using the H53a line parameters given by (Sewilo et al. 2008) we derive 
an LTE electron temperature of T* ~ 7, 600 K, similar to the value for BN and in this case 
consistent with the optically-thin nature of G28.20-0.04N. 

The non detection of H53q: emission from radio source I is consistent with its expected 
large optical depth. The formulation above implies o; ~ 5, and an attenuation factor of 
2/9. This confirms the notion that BN and radio source I arc two sources intrinsically very 
different in nature. This difference is also evident in the brightness temperature of both 
sources. At 7 mm, the brightness temperature of a source is 

(13) 

Using the values of Table 1, we get Tb — 7,800 K for BN, confirming its nature as 
photoionized gas. However, for the radio source I we get Tb — 2,600 K. So, even when 
source 1 seems to be optically thick, its brightness temperature is substantially lower than 
that expected for a photoionized region. Rcid et al. (2007) have discussed as possible 
explanations for this low brightness temperature H~ free-free opacity or a photoionized 
disk. 

Following the discussion of Reid et al. (2007), we consider it unlikely that dust emission 
could be a dominant contributor to the 7 mm emission of BN or Orion 1. A dense, warm, 
dusty disk would be expected to show many molecular lines at millimeter/submillimeter 
wavelengths. While Beuther et al. (2006) and Priedel & Snyder(2008) find numerous, strong, 
molecular lines toward the nearby "hot core", they find no strong lines toward the position 
of Orion I (with the exception of the strong SiO masers slightly offset from Orion I) or BN. 
Also, the brightness temperatures derived by us at 7 mm (7,800 K for BN and 2,600 K for 
source I) are high enough to sublimate dust and suggest that free-free emission from ionized 
gas dominates the continuum emission. Finally, the continuum spectra of BN and of source 
I measured by Plambeck et al.(1995) and Beuther et al. (2006), respectively, suggest that 
the dust emission becomes dominant only above ~300 GHz. 

In the case of source n, no detection was expected given its weakness even in the con- 
tinuum. 



Tb 
K 



0.96 



mJy 



arcsec^ 
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4.3. Spatial Distribution of the H53q; Line Emission 

The H53a line emission in the individual velocity channels shows evidence of structure 
but unfortunately the signal-to-noise ratio is not large enough to reach reliable conclusions 
from the analysis of these individual channels. However, an image with good signal-to-noise 
ratio can be obtained averaging over the velocity range of -21.2 to -1-66.1 km s~^, using the 
task MOMNT in AIPS. This hne image is compared in Figure 3 with a continuum image 
made from the line-free channels. The larger apparent size of the continuum image is simply 
the result of its much better signal-to-noisc ratio. For the total line emission we obtain an 
upper limit of 0"12 for its size, that is consistent with the size of the continuum emission 
given in Table 1. We also show images of the blueshifted (-21.2 to +22.5 km s'^) and 
redshifted (+22.5 to 66.1 km s~^) line emission in Figure 3. The cross in the figure indicates 
the centroid of the total line emission. The centroid of the line emission does not appear to 
coincide with the centroid of the continuum emission and we attribute this to opacity effects. 

An interesting conclusion comes from comparing the total line emission, with the blueshifted 
and redshifted components. The blueshifted emission seems slightly shifted to the SW, while 
the redshifted emission seems slightly shifted to the NE, suggesting a velocity gradient. This 
result supports the suggestion of Jiang et al. (2005) of the presence of an outflow in BN along 
a position angle of 36°. Given the modest signal-to- noise ratio of the data, it is difficult to 
estimate the magnitude of the velocity shift and we crudely assume it is of order one channel 
(~10 km s~^), since most of the line emission is concentrated in the central two channels of 
the spectrum (see Figure 2). The position shift between the blueshifted and the redshifted 
emissions is 0''028±0!'007 (12 + 3 AU at the distance of 414 pc given by Menten et al. 2007), 
significant to the A-a level. Unfortunately, the data of Jiang et al. (2005) does not include 
line observations and there is no kinematic information in their paper to compare with our 
results. 

The small velocity gradient observed by us in BN is consistent with a slow bipolar 
outflow but also with Keplerian rotation around a central mass of only 0.2 Mq. 

5. Conclusions 

We presented observations of the H53a recombination line and adjacent continuum 
toward the Orion BN/KL region. In the continuum we detect the BN object, the radio 
source I (GMR I) and the radio counterpart of the infrared source n (Orion-n) and discuss 
its parameters. In the H53a line we only detect the BN object, the first time that radio 
recombination lines have been detected from this source. The LSR radial velocity of BN 
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from the H53q; line, vlsr — 20.1 ± 2.1 km s~^, is consistent with that found from previous 
studies in near-infared Unes, vlsr = 21 km s~^. We discuss the hne-to-continuum ratio from 
BN and present evidence for a possible velocity gradient across this source. 

LFR and LAZ acknowledge the support of CONACyT, Mexico and DGAPA, UNAM. 
Facilities: VLA 
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Table 1: Parameters of the 7 mm Continuum Sources in the Orion BN/KL Region 





Position" 


Total Flux 






Source 


a(J2000) 


5(J2000) 


Density (mJy) 


Deconvolved Angular 


Size^ 


BN 


05 35 14.110 


-05 22 22.73 


28.6±0.6 


0^'07 ± O^'Ol X 0^'05 ± O^'Ol; 


+ 45° ± 29° 


n 


05 35 14.359 


-05 22 32.78 


3.0±0.6 


< 0^'2 




I 


05 35 14.514 


-05 22 30.57 


14.5±0.7 


0^'09 ± 0('01 X 0('06 ± 0^'02; 


+ 134° ± 44° 



"Units of right ascension are hours, minutes, and seconds and units of declination are degrees, arcminutes, 
and arcscconds. Positional accuracy is estimated to be 0''01. 
''Major axis x minor axis; position angle of major axis. 
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Table 2: Parameters of the 3.6 cm Continuum Sources in the Orion BN/KL Region 





Position" 


Total Flux 




Spectral 


Source 


a(J2000) 


5(J2000) 


Density (mJy) 


Deconvolved Angular Size'' 


Index 


BN 


05 35 14.110 


-05 22 22.74 


4.8±0.1 


0^'17 ± O^'Ol X O^'OS ± 0{'02; + 63° ± 7° 


l.libO.l 


n 


05 35 14.355 


-05 22 32.78 


2.2±0.2 


0^'50 ± 0^'04x < 0{'09; + 20° ± 3° 


0.2±0.1 


I 


05 35 14.514 


-05 22 30.56 


1.2±0.1 


0^'19 ± 0{'04x < 0{'15; + 136° ± 23° 


1.5±0.1 



"Units of right ascension are hours, minutes, and seconds and units of declination are degrees, arcminutes, 
and arcscconds. Positional accuracy is estimated to be 0''01. 
''Major axis x minor axis; position angle of major axis. 
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Table 3: Parameters of H53q; Recombination Line from the BN Object 



Peak Flux 


Half Maximum 


LSR Radial 


Density (mJy) 


Line Width (km s~^) 


Velocity (km s~^) 


10.4±1.1 


39.0±4.9 


20.1±2.1 
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Fig. 1. — Contour image of the 7 mm continuum emission from the BN/KL region. The 
three sources detected are marked with their names. Contours are -5, 5, 10, 20, and 40 times 
0.35 mJy beam^^, the rms noise of the image. The half power contour of the synthesized 
beam (0''17 x 0''15 with a position angle of +17°) is shown in the bottom left corner. 
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Fig. 2. — Spectrum of the H53q; line emission from BN. The dashed hne is the least squares 
fit to the data, whose parameters are given in Table 3. 
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Fig. 3. — Contour images of the 7 mm continuum emission (top left), the average of the 
total (-21.2 to +66.1 km s~^) H53a line emission (top right), the average of the blueshifted 
(-21.2 to -1-22.5 km s~^) HSSa line emission (bottom left), and the average of the redshifted 
(+22.5 to 66.1 km s~^) H53q: hne emission (bottom right) from the BN object. Contours are 
1, 2, 4, 8, and 16 times 1.0 mJy beam~^ for the continuum image. Contours are 1, 2, and 4 
times 12.0 mJy beam~^ km s~^ for the total line image and 16.0 mJy beam~^ km s~^ for the 
blueshifted and redshifted hne images. The cross marks the peak position of the average of 
the total H53q; line emission. The half power contour of the synthesized beam (0'.'17 x O'.'IS 
with a position angle of +17°) is shown in the bottom left corner of the images. 



